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Six-loop massive scheme renormalization group functions of a d = 3-dimensional cubic model 
(J. M. Carmona, A. Pelissetto, and E. Vicari, Phys. Rev. B 61, 15136 (2000)) are reconsidered 
by means of the pseudo-e expansion. The marginal order parameter components number N c = 
2.862 ± 0.005 as well as critical exponents of the cubic model are obtained. Our estimate N c < 3 
leads in particular to the conclusion that all ferromagnetic cubic crystals with three easy axis should 
undergo a first order phase transition. 
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I. INTRODUCTION 



Progress in the qualitative understanding and the quantitative description of critical phenomena to a great extent 
was achieved by the ideas of renormalization group (RG) theory |Q. Only global features of a many-body system 
such as the range of interparticle forces, the space dimensionality d as well as the dimension N and the symmetry 
\Q ' of an order parameter were suggested to be responsible for long-distance and abrupt behaviour of matter in the 
t-H | critical region. As a final step the role of the relevant parameters in microscopic Hamiltonians of various nature 
was represented adequately by effective Hamiltonians used in field theories. While already a vector field theory with 
an isotropic rotationally symmetrical order parameter allowed unified and correct description of a large spectrum 
of critical phenomena, an extension of theories is of special interest since in real substances anisotropies are always 
present Q . For instance, in cubic crystals one expects the spin interaction to react to the lattice structure (crystalline 
anisotropy), suggesting additional terms in the Hamiltonian, invariant under the cubic group. The anisotropy breaks 
' rotational symmetry of the Heisenberg-like ferromagnet and makes the order parameter to point either along edges 
or along diagonals of a cube. The corresponding field theory is defined by a Landau - Ginzburg - Wilson (LGW) 
I ■ Hamiltonian with two (f> A terms of O(N) and cubic symmetry and can exhibit a second order phase transition 
characterized either by spherical or cubic critical exponents. Varying the number of components of the order parameter 
C , N a new crossover phenomenon between these two scenarios takes place at the marginal value N c . 

In the framework of RG theory the critical point corresponds to the stable fixed point of the RG transformation [|l| . 
In a model with competing fixed points the study of domains of their attraction as well as the crossover phenomenon 
is a fundamental problem for universality comprehension. Apart from the academic interest the determination of N c 
can lead to decisive conclusions about phase transition order in a certain class of cubic crystals. For instance, d = 3 
cubic crystals with three easy axis should undergo either a second or a weak first-order phase transition provided N c 
is greater or less than 3 ^ . This argumentation states that the existence of the stable fixed point of a field theory is 
a necessary but not a sufficient condition for a model to exhibit a second order phase transition. If the parameters 
of an microscopic Hamiltonian are mapped in the plain of the LGW Hamiltonian couplings to a point which lies 
outside the domain of attraction of the stable fixed point, the Hamiltonians will flow away to infinitely large values 
of couplings. Such a behaviour might serve as evidence of a weak first order phase transition and this is confirmed in 
some experiments (see Ref. || and references therein). If iV c < 3 for a d — 3 ferromagnet, the new cubic fixed point 
governing the critical regime appears to be inaccessible from the initial values of couplings which correspond to the 
ferromagnetic ordering with three easy axis. It appears that N c is very close to N = 3 and the critical exponents in 
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both regimes are indistinguishable experimentally. In order to calculate the value of N c within field theory one has 
to treat a complex model of two couplings. This is different from a Heisenberg-like ./V-component ferromagnet with 
weak quenched disorder, where the Harris criterion answers the question about the type of critical behaviour pL 

The description of the crossover and the precise determination of its numerical characteristics has been a challenge 
for many RG studies of the anisotropic cubic model. High orders of perturbation theory were obtained for this model 
in successive approximations either in e-expansion with dimensional regularization in the minimal subtraction (MS) 
scheme B or within the massive d = 3 scheme ||. The expressions are available now in the five-loop |Q and in the 
six-loop approximations H respectively. However, divergent properties of the series did not allow their straightforward 
analysis and called for the application of various resummation procedures. For instance, iV c calculated in five-loop 
e-expansion yielded depending on the resummation procedure: N c — 2.958J7), N c = 2.855 ||] and iV c = 2.87(5) [|||. 
Alternative approaches on the basis of the e-expansion lead to N c = 2.97(6) g and to N c = 2.950 On the other 
hand the massive d — 3 scheme RG functions extended for arbitrary N to four loops jnj, yielded iV c — 2.89(2) |ll]] 
(see H for a recent extended review of theoretical determination of N c ). These results suggest that the most reliable 
theoretical estimate is N c < 3. 

However, recent MC simulations [|l2| questioned the values for N c obtained so far. There, considering the finite 
size corrections of a cubic invariant perturbation term at the critical 0(iV)-symmetric point, the eigenvalues uji of 
the stability matrix were extracted. From the estimate u>2 = 0.0007(29) a value of iV c = 3 was concluded. This 
disagreement as well as the crucial influence of the value of N c on the order of the phase transition makes an 
implementation of an alternative method for calculation of N c to be hardly overrated. Recently the massive d = 3 
scheme RG functions of the cubic model were extended to five-loop order [|l3| and very recently the six-loop series @] 
were obtained. The traditional analysis of these series, including an information on large order behaviour of the RG 
functions p), yielded N c = 2.89(4). However let us note that the most accurate estimates of the critical exponents 
of a d = 3 O(N) -symmetric <fi 4 model in a massive field theoretical RG scheme are based on a pseudo-e expansion 
technique Jl4],[l5|]. Up to our knowledge the last method has never been applied to the cubic model jl^]. Therefore 
the main aim of the present paper is to apply the pseudo-e expansion to the up-to-date most precise massive scheme 
RG function || of the cubic model. 

The set-up of the article is as follows. After a brief consideration of the model and renormalization procedure, we 
present the pseudo-e expansion for N c and discuss its properties. Applying Pade- and Pade-Borel analysis we obtain 
precise estimates of N c and compare the result with the corresponding e-expansion. Finally we evaluate the critical 
exponents of a d = 3 cubic system belonging to the new universality class for different values of N > N c and discuss 
the weakly diluted Ising model case N = 0. 



II. PSEUDO-e SERIES AND NUMERICAL RESULTS 



We start from a d — 3 effective LGW Hamiltonian with two couplings at terms of spherical and cubic symmetry: 

r i N / n \ 2 n 

J a=l ' \q=1 / ' a = l 

where ip a (R) are components of a bare iV-component vector field; uo > 0, vq are bare couplings, uiq is a squared bare 
mass being a linear function of temperature. The vicinity of a critical point corresponds to a long-distance behaviour 
of the model (]l|), while ultraviolet divergences of the theory are dealt with by means of an appropriate renormalization 
procedure jj] . In particular, the renormalization of the bare couplings leads to /3 U (u, v) and (3 v (u 7 v) - the so-called (3- 
functions; a renormalization of the bare field and square-field insertion produces y<j,(u, v) and 7^2 (u, v) - the so-called 
7-functions. All these functions depend on renormalized couplings u and v |l8f| . The critical behaviour of the model is 
determined by the infrared stable fixed point u*,v*. It is given by the condition that both /3-functions are zero and all 
the real parts of the stability matrix eigenvalues are positive. The pair correlation function critical exponent 77 and the 
correlation length critical exponent v are obtained via the relations 77 = 7^(u*, v*),l/v = 2 — 7^2 (u*, v*) — 70 (u* , v*). 
The correction-to-scaling exponent to is given by the largest stability matrix eigenvalue in the stable fixed point. 

In the present study we reconsider RG functions of the model ([T]) as they are obtained within massive fixed d = 3 
scheme Q: 
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where (3^ LA ^ 



■ denote the three-loop contributions obtained in Ref. 



|20f|, the four-loop, the recent five-loop 
and the very recent six-loop contributions obtained in Refs. |y]], |l3| and ^respectively. Furthermore, the large- 
order behaviour was established for the RG functions (p|)-(^|) which allowed to apply the refine resummation procedure 
based on Borel transformation with conformal mapping B, . In this way a convergent sequence of approximations for 
N c as well as critical exponents within the cubic universality class were obtained H . 



One possible way of analysis of the massive RG functions 
(resummed) /3-functions (||), (||) 



-(|J) consists in solution of a system of equations for the 



&(«*,«*) =0, 
(3 v (u\v*) =0 



(6) 



to get numerical values of a stable fixed point coordinates u*,v*. Then these numerical values are substituted into 
(resummed) series for 7-functions (||), (||) which lead to the numerical values for critical exponents. As a result the 
final errors for the critical exponents are the sum of the errors of the series for exponents and of the errors coming from 
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u* , v* . To avoid such errors accumulation it is standard now in the analysis of field theories with one coupling [JL4 
to use a pseudo-e expansion fl6fl . Here, we will apply the pseudo-e expansion for a cubic model ([l]) Jlj. The 
procedure is defined in the following way. Let us introduce the functions: 



(3 u (u, v, t) = -u(t -u - -v + ...), 
/3 v (u,v,t) = -v(t - 12 



8 + N 



(7) 



where the "pseudo-e" auxiliary parameter r has been introduced into the /3-functions (Q), (||) instead of the zeroth 
order term. Obviously, (3 u (u,v) = (3 u {u,v,t = 1), /3 v (u,v) = (3 v (u,v,t = 1). Then a fixed point coordinates are 
obtained as series in r. The series for the stable fixed point coordinates u*(t), v*(t) are then substituted into series 
@> @ f° r 7~f unc tfons leading to the pseudo-e expansion for critical exponents. In the resulting series the expansion 
parameter r collects contributions from the loop integrals of the same order coming from both the series of 0— and 7- 
functions. Finally, one puts r = 1. In such a way one gets a self-consistent perturbation theory and avoids cumulation 
of errors originating from different steps of calculation. 

With the above described method we obtain the marginal value N c and the critical exponents in the cubic universality 
class for the model (111). The series for 7V C reads: 



N c = 4 - 4/3r + 0.29042005r 2 - 0.18967704r 3 + 0.19951035r 4 - 0.22465150t 5 . 
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One notes that at least up to the presented number of loops the series does not behave like an asymptotic one with 
factorial growth of coefficients. This can be seen by considering a Pade-table (^|) for N c series (||): 
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Here, a result of an approximant [M/N] is represented as an element of a matrix with usual notation. The approxi- 
mants [0/4] and [1/2] have poles at values of r of the order 1 (at points n = 3.7 and 72 = 1.1 respectively) and thus 
the estimates of N c on their basis are considered as unreliable (they are noted in the table by small numbers). The 
values in the first column of the table are merely the sums of the corresponding number of terms in the expansion (j^) 
and do not diverge. However, the most prominent property of the table is the perfect convergence of the values within 
main diagonals. In particular, the six-loop result of the [3/2] and the [2/3] approximants and the five-loop result of 
[2/21 approximant coincide within the 4th digit and lead to an estimate N c = 2.8616. Though the next order terms 
in (g) could spoil such convergence, it is worth to compare the pseudo-e expansion (||) for N c with corresponding 
e-expansion series which is one loop order shorter Q: 

N c = A-2e + 2.58847559e 2 - 5.87431189e 3 + 16.82703902e 4 . (10) 



The obvious worse convergence properties of the series ( J10| ) lead to a corresponding bad convergence of the values 
in the Pade-table (obtained in Ref. [M). In particular, mere summation of several first terms leads now to diverging 
result: 
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Here, the approximants [0/2], [0/4] have poles close to r = 1 (at n = 2.3, r 2 = 0.9) respectively and thus are 
unreliable. 

In order to take into consideration possible factorial divergence of the pseudo-e expansion (||), as a next step we 
apply to the series (^) the Pade-Borel resummation procedure. The Pade-Borel resummation of the initial sum 
N c (t) consists of the following steps: i) construction of Borel-Leroy image of N c (t); ii) its extrapolation by a rational 
Pade-approximant [M/N}(rt) and iii) definition of a resumed N^ es (r) by the integration J^ dtexp(—t)t p {M/N](Tt), 
where p is an arbitrary parameter entering the Borel-Leroy image [pTJ . One possibility to fix p is to require fastest 
convergence of the resulting values, given by the diagonal approximant resummation similar to Pade-analysis (|^). 
However the convergence of these values appears to be almost independent of p. On the other hand approximants 
possessing poles on the positive real axis are considered as unreliable and we can equally well choose p to provide a 
minimal number of such divergent approximants. For instance, for p > 4 the imaginary part is smaller then 10 -10 
in the "bad" [1/4] and [3/2] approximants and therefore can be neglected. Processing the series (||) for p = 4 as 
described we obtain the results presented in (|l2|). Here, one encounters only one unreliable approximant [1/2] which 
is again denoted with small numbers. This analysis yields N c = 2.862±0.005, where error bar stems from the maximal 
deviation between the six and the five-loop results for arbitrary p between and 10. 
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It is obvious that other values of interest such as fixed point coordinates and critical exponents can also be obtained 
within the pseudo-e expansions. For instance, for different values of N we obtain the following expressions for the 
critical exponents 7 of the susceptibility, v of the correlation length, and uo of the correction-to-scaling: 
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One can determine the cubic model critical exponents of the new universality class on the basis of the expansions jl^ ) 
in the same manner as for the expansion (|^) of N c . However, the expansions for the combination 1/7 and 1/v — 1 
appear to have better convergence properties and all values are obtained on their basis. The Pade- and Pade-Borel 
analysis lead to the critical exponents values as they are given in the table [j] in the last column (here, we do not show 
intermediate results similar to (|9])-(|l2"|)). The error bars for the critical exponents given in the Table Q were obtained 
from the maximal deviation between the six- and the five-loop results among all deviations for the parameter p value 
< p < 10. The error bars within the pseudo-e expansion are typically much smaller than those based on other 
methods. The reason is that (3- and 7-functions contribute in a self-consistent way into the pseudo-e expansion series 
for critical exponents. 

To compare our results we represent in the table | the values obtained from the five-loop e-expansion by means 
of a modified Borel summation |E| and of the Borel transformation with conformal mapping Q (the corresponding 
citation in the table is primed). Critical exponents from the four-loop fixed massive d = 3 scheme with an application 
to the RG functions of Pade-Borel resummation jll|] and the results from the six-loop RG functions resummed by Borel 
transformation with conformal mapping ||] are given in the table. Recently, the modified Pade-Borel resummation |22] ] 
to the six-loop RG functions l| of the cubic model. These data are also displayed in the table. 

The error bars for the values of the Refs. || and |24| were obtained from the condition of the result stability in successive 
approximation orders. However, the numerical value of the fixed point coordinates were substituted into the expansions 
for the 7-functions (||), (||). To this end the most reliable numerical values of the stable fixed point coordinates were 
substituted into the resummed 7-functions and then an optimal value of the fit parameter in modified Pade-Borel 
resummation p4| | and two fit parameters in the conformal mapping procedure j|] were chosen. The deviations between 
five- and six-loop results obtained within the resummation procedure with optimal fit parameter(s) value gave the 
error interval. To complete the list, we show the value of uo for N = 3 obtained in Ref. on the basis of Borel 
transformation involving knowledge on RG functions large-order behaviour. Let us note that for finite values of N 
our data for 7 and v interpolate the results of the minimal subtraction ( |22| , || ) and the massive scheme ( |llj , [g ) , 
though the values are closer to the last. On the other hand, our method gives smaller values for uj in comparison with 



has been applied |24 
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other methods. We note as well that passing from the four-loop to the six-loop |24j approximation in frames of 
a massive d — 3 approach shifts the numerical values of critical exponents towards our data. 

In the limit N — > oo the critical properties of the cubic model reconstitute those of the annealed diluted Ising model pq] 
where Fisher renormalization for critical exponents holds [ p6| . In particular, based on the recent RG estimates for 
the critical exponents of the pure Ising model a = 0.109 ± 0.004, v = 0.6304 ± 0.0013, 7 = 1.2397 ± 0.0013 Q, one 
obtains the values v — 0.708 ± 0.005, 7 = 1.391 ± 0.008 for the Ising model model with annealed disorder. The last 
values agree very well with our results of the last row of the table |. Moreover, they are in very good agreement with 
other data of the table. 

It is worth to note here that the RG series for the cubic model allow to reconstitute the functions which describe the 
Ising model with the other type of randomness. By substitution N = one reconstitutes the weakly diluted quenched 
Ising model (RIM) j2^]. In this case, however, the pseudo-e expansion in t degenerates into a -^/r-expansion for 
the same reasons as the e -expansion for the RG functions degenerates into a y^-expansion p7| . Moreover, our 
calculation show that an expansion in \fr is numerically useless as this was shown for y^-expansion [^|, [Io|, S . 
This can be regarded as an evidence of the Borel non-summability of the RIM RG functions. Since the asymptotic 
properties of the series still are not proven despite of noticing their divergent character |^,|3(|, the RG functions of 
RIM as series in renormalized couplings used to be treated by means of Pade-Borel or Chisholm-Borel resummations 
(see e.g. Jl0| , j3l|). The first of the stated methods was recently applied to study the five- loop RG functions 
of RIM ]13| . The analysis allowed the authors to obtain the five-loop estimates for the RIM critical exponents. 
Extending the analysis of Ref. |l3| to the six-loop order reveals the wide gap between five- and six-loop fixed point 
coordinates. This leads to an inconsistency of the six-loop values of critical exponents compared with the five-loop 
results of Ref. |l3|| . However, the analytical solution of a toy d — RIM showed its free energy to be Borel summable 
provided that resummation is done asymmetrically: resumming first the series in the coupling u and subsequently the 
series in v |50| . The corresponding resummation applied to d = 3 RIM massive scheme RG functions allowed precise 
determination of the critical exponents |52| . 



III. CONCLUSIONS 



In the present paper we studied the critical properties of a cubic model associated with cj) 4 -teims of spheric and cubic 
symmetry of the LGW Hamiltonian. In particular, we were interested in the crossover between 0(7V)-symmetric and 
cubic behaviour which occurs at a certain value N c of order parameter components number. Recently, five- |l3| and 
six-loop Q order RG functions were obtained for the cubic model within massive d = 3 scheme ||. We applied the 
pseudo-e expansion to their analysis. This method is known as a standard one for the (3(iV)-symmetric model 
analysis and leads to the most accurate values of critical exponents [Il5l. Here, to our knowledge, it has been applied 



to the cubic model for the first time 17 



The pseudo-e expansion for N c appears to have much better convergence properties then the corresponding e- 
expansion (c.f. Pade-tables ^ and (0)). This provides very good convergence of its Pade-analysis ([|). The last 
together with the refined Pade-Borel analysis yields the best estimate N c — 2.862 ±0.005 of the paper. Our conclusion 
N c < 3 means in particular, that all ferromagnetic cubic crystals with three easy axis should undergo a first order 
phase transition || . 

We obtained the values of cubic model critical exponents in the new universality class in pseudo-£ expansions with 
the results given in table [|. In the N — > 00 limit our data reproduce the critical behaviour of an annealed weakly 
diluted Ising model |25). The N — > limit, corresponding to a quenched weakly diluted Ising model [^Tj, however does 
not yield reliable results in pseudo— y/e expansion. Within a traditional d = 3 massive technique the resummation 
of the RG functions by means of the convenient Pade-Borel analysis reveals a gap between five- and six-loop fixed 
point coordinates. This leads to an inconsistency of the obtained critical exponents values compared to the declared 
in Ref. |l3[ |. Let us note, however, that recently reliable values have been obtained |32| by a resummation method 
which treats the couplings of the RIM model asymmetrically 130] . 
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TABLE I. Our data for the critical exponents of the cubic model (last column) in comparison with other results. See the 
text for a full description. 



N 



22 



[HI 



[24\ 



this study 





7 


1.3746 ± 0.0020 


1.377(6) 


1.3775 


1.3850 ± 0.0050 


1.390(12) 


1.387 ±0.001 


3 


V 


0.6997 ±0.0024 


0.701(4) 


0.6996 


0.7040 ± 0.0040 


0.706(6) 


0.705 ± 0.001 




UJ 


0.8061 pi 


0.799(14) 


0.7786 


0.7833 ± 0.0054 


0.781(4) 


0.777 ± 0.009 




7 


1.4208 ± 0.0030 


1.419(6) 


1.4028 


1.4074 ± 0.0030 


1.405(10) 


1.416 ± 0.004 


4 


V 


0.7225 ± 0.0022 


0.723(4) 


0.7131 


0.7150 ±0.0050 


0.714(8) 


0.719 ±0.002 




UJ 




0.790(8) 




0.7887 ± 0.0090 


0.781(44) 


0.777 ± 0.002 




7 


1.4305 ± 0.0040 




1.4076 






1.417 ±0.006 


5 


V 


0.7290 ± 0.0016 




0.7154 






0.720 ± 0.004 




UJ 












0.773 ± 0.003 




7 


1.4322 ± 0.0040 




1.4082 






1.417 ±0.009 


6 


V 


0.7301 ± 0.0016 




0.7157 






0.718 ± 0.003 




UJ 












0.771 ± 0.005 




7 




1.422(6) 


1.4074 


1.4068 ± 0.0030 


1.404(10) 


1.407 ± 0.008 


8 


V 




0.723(2) 


0.7153 


0.7143 ± 0.0035 


0.712(6) 


0.715 ±0.003 




UJ 




0.786(6) 




0.7955 ±0.0150 


0.775(88) 


0.770 ± 0.009 




7 


1.3993 ± 0.0020 


1.399(8) 




1.3962 ± 0.0040 


1.396(14) 


1.395 ± 0.006 


oo 


V 


0.7108 ±0.0010 


0.711(2) 




0.7094 ± 0.0030 


0.708(8) 


0.708 ± 0.001 




UJ 




0.802(18) 




0.7986 ± 0.0200 


0.790(18) 


0.775 ± 0.020 
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